GAUGE APPROACH TO THE SPECIFIC HEAT 
IN THE NORMAL STATE OF CUPRATES 
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Many experimental features of the electronic specific heat and entropy of high Tc cuprates in the 
normal state, including the nontrivial temperature dependence of the specific heat coefficient 7 and 
negative intercept of the extrapolated entropy to T = for underdoped cuprates, are reproduced 
using the spin-charge gauge approach to the t — J model. The entropy turns out to be basically 
due to fermionic excitations, but with a temperature dependence of the specific heat coefficient 
controlled by fluctuations of a gauge field coupling them to gapful bosonic excitations. In particular 
the negative intercept of the extrapolated entropy at T = in the pseudogap "phase" is attributed to 
the scalar component of the gauge field, which implements the local no-double occupancy constraint. 

PACS numbers: 71.10.Hf, 11.15.-q, 71.27.-|-a 



The low-temperature electronic entropy of high Tc 
hole-doped cuprates in the normal ( "metallic" ) state ex- 
hibits a behavior rather unusual for a metal: the specific 
heat coefficient 7 which should be constant shows a non 
trivial temperature dependence [l| and even more spec- 
tacularly the T = intercept of the entropy, extrapolated 
from the approximate linear behavior at moderate tem- 
peratures, turns out to be negative in the underdoped 
region [2|. In this paper we apply the spin-charge gauge 
approach developed in [3|, Ij, y, l6| to extract from the 
two-dimensional t-J model the low temperature entropy 
and specific heat in the normal state and compare our 
results with the experimental data, in particular show- 
ing how this approach can explain the peculiar behavior 
mentioned above. 

Let us first outline the main features of the experi- 
ments [l|, y, lZI , following Ref. [l|. In the normal state of 
La2_5Sr5Cu04 (LSCO) and YBazCuaOg-Hx (YBCO) the 
electronic specific heat coefficient 7 = C^^ jT as a func- 
tion doping concentration 8 and temperature up to 250- 
300 K exhibits the following behavior. For strongly un- 
derdoped samples 7 ((5, T) -- A7^^-HB((5, T) where A7^^ 
is a Low- Temperature upturn, B(S,T) is slowly increas- 
ing roughly linearly in S and T with almost (5-independent 
slope. At higher dopings the Low- Temperature upturn 
disappears and the samples become superconducting. 
The increasing part remains, but it saturates to a broad 
maximum at T* followed by a slow decrease in T. T* 
roughly coincides with the pseudogap temperature T* 
identified by the inflection point in the in-plane resistiv- 
ity, as can be checked using the data on curvature of re- 
sistivity [8|, see Fig 1. In the decreasing region 7 becomes 
almost (5-independent. Similar features are exhibited, in 
the appropriate doping range, also by Bi2Sr2CaCu208+i^ 
(Bi2212) data [9|], so they can thus be considered as 
rather generic in cuprates and therefore it is reasonable to 




FIG. 1: The dots mark the inflection point on resistivity (T*) 
as given in Q] on 7 data in YBa2Cu306+:E taken from [y]; 
semiempirically S ~ 0.2a;, see [2|. 



explore a physical interpretation in terms of the physics 
of doped Mott insulator described by the t — J model. 

An attempt to discuss the specific heat coefficient in 
terms of SU{2) slave boson theory [10[ of the t — J model 
appears in Ref I 111, obtaining for the gauge contribution 
7 « — TlnT and some data are fitted with this formula. 
Further attempts to understand the behavior of the spe- 
cific heat can be found in [l2] [l3[. 

In spite of the complex features of "f[5,T) described 
above, the entropy S{5, T) of the normal state exhibits 
a simpler continuity: for low dopings it is approximately 
linear in T with increasing slope as 5 increases, at higher 
T the slope becomes almost i5-independent with increas- 
ing but negative intercept at T = 0. At higher dop- 
ings {5 > 0.19) also the intercept becomes almost S- 
independent and approximately 0. On the basis of the 
above continuity it was argued in Ref. [l] that in YBCO 
the spin excitations may be relevant for all S and the 
above results are better described as a modification of 
the low-energy spin spectrum as 5 changes than by a 



simple band model. The spin-charge gauge approach ap- 
pears to partially substantiate such claims. The negative 
intercept of the entropy suggests a negative contribution 
to entropy of a "constraint-field" , which acts reducing 
the low-energy degrees of freedom or more precisely re- 
moving them from the temperature/energy region con- 
sidered. This was proposed in the analysis of the ther- 
modynamics of the t — J model performed in 1J| within 
the slave-boson approach. In fact, a negative contribu- 
tion to entropy naturally arises in a gauge approach from 
the scalar component of the gauge field (in the Coulomb 
gauge) enforcing "Gauss law" . Let us now sketch the ba- 
sis of the spin-charge gauge approach and its application 
to the computation of entropy and specific heat. 

This approach assumes as a (simplified) model for CuO 
layers in high Tc cuprates the 2D t-J model with t/J^3. 
Neglecting t', t", details of Fermi surface (FS) are lost but 
the analysis is simplified, hopefully retaining the basic 
relevant features. The model is treated in an "improved 
Mean Field Approximation" (MFA) via a gauge theory 
of spin-charge decomposition, obtained by gauging the 
global spin and charge symmetries of the model [l5| . This 
gauging is obtained introducing spin and charge Chern- 
Simons gauge fields. The nice feature of introducing these 
gauge fields is the possibility of a more flexible treat- 
ment of charge and spin responses within a a spin-charge 
decomposition scheme. In the end they will disappear 
from the game in MFA, but leaving behind sign of their 
presence crucial for the low-energy physics, as discussed 
below. 

The basic fields adopted in this approach for the spin- 
charge decomposition of the t — J model are a charged 
spinless fermion, the holon [16|, a neutral spin 1/2 bo- 
son of a non-linear a (CP^) model, the spinon, and a 
slave-particle gauge field (not to be confused with spin 
and charge Chern-Simons gauge fields). The spin-gauge 
field in MFA attaches spin vortices to the empty-site po- 
sitions. The spinous moving across this gas of vortices 
acquire a mass gap, with a theoretically derived doping 
dependence, rus ~ ■\/|(51n^ consistent with AF correla- 
tion length at small 5 derived from neutron experiments 
[l7|. In MFA at low temperature and small doping con- 
centration, in the parameter region to be compared with 
the "pseudogap phase" (PG) of the cuprates, the holons 
move in a statistical magnetic field with flux tt per pla- 
quette generated by the charge-gauge field. This "phase" 
shares some similarity with the 7r-flux phase appearing in 
the slave boson formalism [l8|. 

Around the pseudogap temperature T* the 7r-flux lat- 
tice "melts" and we enter in the "strange metal phase" 
(SM), at higher 5 ov T , see Fig. 2. Notice that, since 
only holons are involved and not full electrons, this is not 
a true phase transition as the one appearing in the DDW 
formalism [l9|. 

In PG, as a consequence of the 7r-fiux, the holons are 
converted via Hofstadter mechanism into two species of 
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FIG. 2: Qualitative phase diagram with in grey the "phases" 
considered in the paper 



Dirac fermions with small Fermi surface {ep ~ t5) cen- 
tered at the four nodes ( ±^, ±^ J, whereas in SM they 

exhibit a large Fermi surface {ep ~ t{\ — <?)), as expected 
from band structure calculations. A direct evidence of 
the small FS in PG might come from recent experiments 
on Shubnikov-de Haas oscillations [20|. 

Holons and spinous are gauge-invariantly coupled by 
a C/(l) slave-particle field, A, whose low energy effective 
action is obtained upon integration of the matter fields. 
As a consequence of the finite FS of holons the transverse 
gauge propagator exhibits a Reizer singularity [2l| which 
dominates at large scales: for small g,a;,w/|g| 



{A^Ai_){u, q) - {-xW + ii^^/W 



(1) 



where A±^ is the transverse component of A, x is the 
diamagnetic susceptibility and k the Landau damping. 
Both x~^, K and the holon mass mh are ~ 5 in PG and 
'^ 1 — (5 in SM. The scalar component Aq has a low energy 
propagator given by 
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where rriQ is a thermal mass generated by the spinous 
and H the Heaviside step function. In view of the con- 
stant term in ^ the interaction mediated by Aq is short 
ranged, hence subleading at large distance w.r.t. the in- 
teraction mediated by A^. However, taking into account 
the renormalization of the transverse contribution dis- 
cussed below, it gives the dominant contribution to 7 in 
PG for T of the order of the holon Fermi temperature, 
which in this "phase" is rather small, of the order of a 
hundred K at low dopings. In the gauge correlator the 
momenta extend up to an UV cutoff A ^ J. 



As discussed in [J, |l8|, |2^ at finite temperature T 
the typical momentum scale of the transverse gauge 
fluctuations is given by the anomalous skin momentum 
Qo = {^T/xy^^- The interaction of holon and spinon 
modifies the gauge propagator, inducing a cutoff of the 
infrared momentum singularity, for \q\ < Qq, replacing 
k|(7|~^ by the sum of the conductivity of the spinon- 
gauge and holon-gauge subsystems, denoted by as and 
ah, respectively [23|. The splitting in high and low mo- 
menta contributions is explicitly realized in the calcula- 
tions with a sharp cutoff at |g| = (Qo, where ^ ~ 0.4 25 1 
and assuming a renormalization of k at high momenta 
accordingly. 

To extract the entropy and the specific heat we start 
by computing the free energy F; then S — —^ and 

7 = ^^7^1 where we only differentiate w.r.t. the explicit 
dependence on T, thus ensuring S{T = 0) = 0. Within 
our approach F is the sum of four terms: the contribu- 
tion of free spinous. Fg, of free holons, Fh and the fully 
renormalized contribution of transverse and scalar gauge 
fluctuations, Fj_,Fq. Since the spinons arc massive, the 
T-dependence of Fg is negligible for T lower than the 
spinon gap, which we estimate of the order of few hun- 
dreds K. For the holons we have the standard result (for 
each holon species) 



Sfi « cruhT 



(3) 



where the phenomenological constant c accounts also for 
the eccentricity of the FS due to neglected t',t" terms. 
A comparison with [26| yields c « 3. We estimate 
Ftt,tt =-L,0 following Ref. [13 if we denote by Dj the 



fully renormalized retarded Green function of the gauge 
held 

duJC0th{u;/2T) [ d\aictaii{^^^^^^^^). (4) 

We remark that in the calculation of Fq in [l^ a so- 
phisticated "ad hoc" regularization was needed because 
the euclidean scalar correlator vanishes in the limit of 
infinite frequency, thus making impossible a direct ap- 
plication of ^-function regularization, since both spinons 
and holons are gapless in the slave-boson approach. This 
problem does not arise here due to the constant term 
in the scalar correlator caused by gapful spinons. The 
key result of [1J| that Fq and F±_ have opposite sign is 
recovered here as a consequence of the opposite relative 
sign of the real and immaginary part of D^ for scalar and 
transverse components, as follows from eqs. ([T]) and ^. 

The dominating contributions [28| to entropy of trans- 
verse gauge fluctuations turn out to be, up to logarithmic 
corrections. 
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where a = ah + a, ^^ n^p + xm^T-^ 0. The T^/s be- 
haviour is the standard one for 2D clean electrodynamics 
27| . The second contribution in SM comes from "small" 
momenta and is negligible in PG. The leading scalar con- 
tribution is given by 



VF 



(6) 
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(5) 



where vp is the holon Fermi velocity. One can verify that 
increasing T first the transverse then the scalar contri- 
bution dominates in PG, whereas in SM, in the tempera- 
ture range considered the transverse contribution always 
dominates. 

/^From equations ([3]), (O and ([6|) one can easily derive 
the following consequences for the interpretation of ex- 
perimental data within the spin-charge gauge approach: 

1) The approximately linear behaviour of S is basically 
due to the holons, although it is renormalized by gauge 
fluctuations. The increase of the slope at low 5 and its 
saturation at higher T or 5 are due to the transition from 
•nih '^ S characteristic of PG to nih '^ 1 — d characteristic 
of SM. 

2) The negative intercept of entropy is due to the scalar 
gauge contribution ([6|) in PG, negative in agreement with 
the general ideas discussed in the introductory section. 

3) The upturn Aj^'^ in 7 (and the analogous more ev- 
ident in S/T [2] ) is due to the contribution of transverse 
gauge fluctuations in PG. Their contribution in SM yields 
the decrease above T* which we identify as the PG-SM 
crossover. These enhancement of entropy are due to the 
presence of the gapless transverse gauge mode. 

4) The approximately linear increase of 7 in (5 and in 
T, with almost S independent slope, is due to the second 
term in ^ , presumably together with a pairing con- 
tribution (see below). This behavior replaces the linear 
slope due to the contribution of AF spin waves in the 
pristine material, now removed by the spinon gap. 

More concretely a comparison between theory and ex- 
periments is summarized in figures 3 and 4, where in 
the inset the experimental data for selected dopings are 
plotted only for the region of parameters discussed above, 
where a comparison is meaningful. The theoretical curves 
have been obtained with the same value of the parame- 
ters used in [J, |5|, l6| and by substituting in a the ex- 
pression for the conductivities derived there. The value 
of vp in PG extracted from the slope of experimental 
7 in the increasing range compared with the theoreti- 
cal expression derived in from ([6]) turns out to be of the 
order of magnitude of the electron Fermi velocity found 
experimentally in ARPES [26| and of the holon Fermi 
velocity used in the calculation of transport properties 
in [J, y, although 2-3 times smaller. Presumably this 
is due to a further contribution of pairing , as in the 
preformed pair (see e.g. 18|, l29|) or fiuctuating phase 
superconductor [30| approaches, which would yield a T- 
incrcasing density of states, not taken into account in the 
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FIG. 3: The calculated temperature dependence for S (above) 
and 7 (below) for selected dopings 5 in both PG and SM 
"phases" , in comparison with experimental data for analogous 
in-plane doping concentrations in YBa2Cu3 06+a; ), taken 
from PI (inset; labels show x). 



present simplified treatment. This increase should also 
account for the smooth transition from PG to SM dis- 
cussed above in item 1). However a depletion mechanism 
of density of states (DOS), alone seems to be unable to 
reproduce the experimental behavior of 7, because lower- 
ing T the curves at different dopings should converge near 
T = with increasing slope as b increases (see e.g. [9[, fig 
12), whereas experimentally they are basically parallel at 
moderate temperatures with an upturn at low T in non- 
superconducting samples. The presence of two distinct 
effects appears consistent also with recent experimental 
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data on specific heat in Bi2Sr2_a;Laa:Cu06 

As one can see, even in our simplified treatment a qual- 
itative agreement of the behavior both in T and 5 is found 
deeply in PG and SM in the region of validity discussed 
above, reproducing the features discussed in the intro- 
ductory section. The derived behavior of the entropy in 
SM is also consistent with the numerical data extracted 
from the t-J model with the Lanczos method [34I or high 



temperature expansion [lj| in the overlapping range of 
temperature. 

Summarizing, the reasonable agreement of doping- 



FIG. 4: Detail: Calculated S in pseudogap "phase", in com- 
parison with experimental data on YBCO (inset) taken from 
[1[ with the same notations of Fig. 3 



temperature dependence of experimental data and the 
theory discussed above suggests that one can interpret 
the low-temperature electronic specific heat in the nor- 
mal state of hole-doped high Tc cuprates basically as due 
to fermionic holons, with small FS in the "pseudogap" 
and large FS in the "strange metal phase", but gauge 
fluctuations determine the variation in the T-dependence 
of the specific heat coefficient . Gauge fluctuations them- 
selves reflect the changes in the spectrum of renormalized 
spin excitations and charge carriers. 

In this approach the rather surprising negative inter- 
cept for the entropy is due to the negative contribution 
of scalar gauge fluctuations. Since the gauge field is just 
a consequence of the Gutzwiller projection, this interpre- 
tation in terms of a "constraint" gauge field suggests that 
this peculiar feature is rooted in the no-double occupancy 
condition, hence a basic characteristic of the low-energy 
description of doped Mott insulators. The phenomenon 
appears in PG due to the "effective" rather small Fermi 
temperature, a consequence of the 7r-flux, characteristic 
of a 2D system, suppressing the positive transverse con- 
tribution. 

Furthermore, in this approach the gauge field is the 
"glue" binding spinon and holon into an electron reso- 
nance [a, ISj . The above interpretation then strengthens 
the idea of a composite nature of the charge carriers in 
the cuprates, an idea which is also strongly suggested by 
the metal-insulator crossover, as discussed in detail in 

n. 
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